ANOMALOUS TAG DIFFUSION IN THE ASYMMETRIC EXCLUSION MODEL 
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Anomalous behavior of correlation functions of tagged particles are studied in generalizations of 
the one dimensional asymmetric exclusion problem. In these generalized models the range of the 
hard-core interactions are changed and the restriction of relative ordering of the particles is partially 
brocken. The models probing these effects are those of biased diffusion of particles having size S = 
0, 1, 2, ... , or an effective negative "size" S = — 1, —2, . . . , in units of lattice space. Our numerical 
simulations show that irrespective of the range of the hard-core potential, as long some relative 
ordering of particles are kept, we find suitable sliding-tag correlation functions whose fluctuations 
growth with time anomalously slow (t 3 ), when compared with the normal diffusive behavior 
These results indicate that the critical behavior of these stochastic models are in the Kardar-Parisi- 
Zhang (KPZ) universality class. Moreover a previous Bethe-ansatz calculation of the dynamical 
critical exponent z, for size S > particles is extended to the case S < and the KPZ result z = | 
is predicted for all values of S £ Z. 
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Dynamic and equilibrium properties of lattice gases, 
i. e., systems involving randomly moving particles with 
hard-core interactions, have received much attention 
within the last decades, mainly due to their intrinsic and 
nontrivial many body behavior. Several important re- 
sults have been obtained for such systems of particles. In 
particular the study of tagged-particle diffusion in hard- 
core lattice gases with nearest-neighbor hopping revealed 
interesting nontrivial many body behavior In one 

dimension the rate of diffusion depends strongly if the dif- 
fusion is symmetric or not, i. e., has or not a prefered 
direction. In the symmetric case, the root-mean-squared 
displacement of a tagged particle increases as 1 1 during a 
time interval t, behavior that is anomalously slow when 
compared with the asymmetric case where this quantity 
exhibits the standard diffusive behavior t?. Majumdar 
and Barma jnj have shown that even in the asymmet- 
ric case it is possible to find anomalous behavior of spe- 
cial correlations of tag particles (sliding-tag correlations) , 
where the difference of the tag coordinates of the parti- 
cles changes in time with a fixed tag velovity vtag- For 



arbitrary values of Vtag, except for v tag 



these 



correlations increase with a typical diffusive behavior t 2 , 
while at the special value vtag = v^ ag , that depends on 

the density of particles, they increase anomalously as 43. 
This last behavior can be understood by mapping the 
asymmetric diffusion model into the particle height in- 
terface model, whose fluctuations in the continuum limit 
are governed by the Kardar, Parisi and Zhang (KPZ) 
model Q. 
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From the above results (T^j it is clear that the relevant 
parameters ensuring anomalous behavior for correlations 
of tag particles are: the hard-core effect (that implies, 
due to the topology of a one dimensional lattice, a fixed 
relative ordering of particles), the right-left asymmetry 
of the diffusion, and the velocity in the tag space where 
the correlations are measured. 

All these results are obtained by imposing that the tag 
particles has a unity size, in terms of lattice units, or 
equivalently the particles have a fixed relative order and 
two particles at the same site have an infinite hard-core 
repulsion. The aim of this paper is to obtain the anoma- 
lous behavior of these sliding-tag correlations in the case 
where we increase or decrease the range of exclusion of 
molecules as well we relax the constraint of relative order- 
ing of the molecules. The variation of the exclusion effect 
will be done by considering the asymmetric diffusion of 
particles with arbitrary integer sizes S = 0, 1, 2, ... , and 
the effect of ordering will be verifyed by considering a 
generalization of the model where two consecutive or- 
dered particles, can violate partially their ordering. This 
last model, as we shall see, corresponds to a generaliza- 
tion to the case where the particles in the asymmetric dif- 
fusion have an effective negative "size" S = —1,-2, . . .. 

We consider initially a one dimensional periodic lat- 
tice of N s sites contaning N p particles of a fixed size S 
(S = 0, 1, 2, ... ). These particles occupy S consecutive 
sites, having a hard-core repulsion of range S. In the case 
5 = the range of infinite repulsion is zero and we may 
put an arbitray number of particles in a given site. The 
dynamics of the asymmetric diffusion is such that a par- 
ticle is selected at random and with probability p (q) the 
particle attempts to move by one unity of lattice spacing 
to the right (left) . The movement is accepted if the hard- 
exclusion allows the final configuration. The unit of time 
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is given by N p attempts of motion. For a given number 
of particles N p the number of holes (positions where the 
particles are allowed to move) is given by N s — SN p and 
the mean velocity v p of the particles is given by 



(p-q)(N s -SN p ) 
N P + N S - SN P 



= (p-q)(l- X ), (1) 



where \ = N p /(N S — (5 — 1)N P ) plays the role of an 
effective density of particles in the system. A given con- 
figuration at time t is given by the lattice coordinates of 
the particles {y(n, £)}, where n = 1, 2, , . . . , N p are the 
sequential labels indexing the particles. Since indepen- 
dently of the size of molecules the asymmetric diffusion 
process does not change their order (even for the size 
5 = particles), we may chose this ordering such that 
y(n + l,t) > y(n, t) + 5, and the periodic boundary con- 
dition translate into y(n ± N p , t) = y(n, t) ± N s . A cor- 
responding interface fluctuation model can be defined by 
considering y(n, t) as the local height of the surface at the 
horizontal cordinate n. In this model, at each time step, 
with probability p (q) y(n,t) increase (or decrease) by 1, 
and the motion is performed if y(n + 1, t) — y{n, t) > 5 
(or y(n,t) — y(n — l,f) > 5). In the steady state, for 
p > q (or p < q) the surface grows (or decreases) in the 
vertical direction, with velocity v p . 

In order to search for anomalous behavior of correla- 
tions we consider a set of correlations among the coordi- 
nates {y(n t ,t)} 7 of the tagged particles with label {n t } 
at time t and the cordinate y{n$, 0) of a given particle hq 
at initial time t = 0. The label difference of the tagged 
particles nt — no = Vtagt changes with time, and for con- 
venience we write the tag velocity as Vtag — pv p b, where 
p = Np/N s is the density of particles and b = b(S,x) is 
a real parameter. The mean values of these coordinates 



< y(nt,t) - y(n ,0) > 



< y(n t , 0) - y(n , > +v p t 

(Vp + V ta g/p)t = V F t (2) 



increases with the frame velocity vp 
sliding-tag correlation functions tpT[- 
are given by 



v p (l + b). The 
! we considered 



cr(vtag,t) 2 =< [y(n t ,t) -y(n ,Q)} 2 > -v 2 F t 2 . 



(3) 



We have measured these correlations by extensive Monte 
Carlo simulations by varying the asymmetry parameters 
p, q, the density of particles p, and the size 5 of the par- 
ticles. We verified that irrespective of the size 5 of the 
particles, as long the diffusion is biased (p ^ q), anoma- 
lous behavior happens for special values of the tag ve- 
locity Vtag- The value of Vtag depends on the density p, 
the asymmetry parameter (p — q) and the size 5 of the 
particles. Examples of these measurements, for the cases 
where 5 — 2 and 5 = are shown respectively in Fig. 1 
and 2, where the time dependence of the correlation (||) 
are shown for several values of the tag velocity. In Fig. 
1 we consider p = 0.75, q = 0.25 and density p = 1/5 in 
a lattice of N s = 150 000 sites while in Fig. 2 p = 0.75, 



q = 0.25 and p = 1/4 in a lattice of N s = 120 000 sites. 
The average slope of the curves on these figures shows 
that apart from crossover effects, due to the finite size 
of the lattice, we have a normal diffusive behavior for 



arbitrary values of Vtag, except at Vtag 



Hag 



where b c = 0.4167 for S = 2 and b c = 0.2 for 5 = 0. 
At those critical values the increase of the correlations 
is anomalously slow and our numerical results indicates 
a ~ t" where v = 0.32 ± 0.01, for the 5 = particles and 
v = 0.33 ± 0.01 for the 5 = 2 particles. 

A constraint on the dependence of the parameter 
b c = b c (S,x) can be derived by generalizing the argu- 
ments presented in [ p"T[ for unity size particles. Instead 
of describing the asymmetric diffusion in terms of the 
size 5 particles we consider the diffusion of holes in the 
reverse order. This is possible since the holes does not 
change their relative order and behave as ordinary size 
5 = 1 particles traveling on an effective lattice with 
N s - N p (S - 1) = 1 - x(5 - 1) sites (allowed posi- 
tions). The mean velocity of holes is Vh — 5x(<Z — p), 
where the factor 5 appears because in their motion the 
holes jump particles with size 5 We can consider 

an equivalent sliding-tag correlation function (^) for the 
holes, with the tag velocity v^ ag = (1 — Sp)b h {\)vh and 
the frame velocity v F = (1 — b h (x))vh- The correspon- 
dence between 6(5, x) an d b h (x), describing the same 
fluctuations in both correlations (particles or holes) are 
obtained by imposing the equality |ll| vp — v F . At 
the critical point hole and particle correlations should 
reveal the same anomalous behavior. Taking into ac- 
count that the holes behave as ordinary 5=1 particles, 
but traveling 5 times faster in an effective lattice with 
N s — (5 — 1)N P sites, we have by comparing the tag ve- 
locities, (A^ s - (5 - l)N p )Sb*(x) = NM1, 1 - x)- As a 
consequence we obtain the general relation 



X (l - x(l - S))b c (l, 1 - x) + (1 - x)b c (S, x) = 
l + (S- l)x, (4) 

that recovers the known relation in the particular case 
where 5 = 1 @. 

Our numerical experiments show that the critical val- 
ues 6 C (5, x) do not depend on the asymmetry parameter 
as long the diffusion is asymmetric (p ^ q). Moreover 
our results are consistent with the conjecture 



b c (S,x) = ( 



S 



i-x 



+ 1 - S) X - 



(5) 



The agreement of the above conjecture was tested ex- 
haustively in numerical simulations for several values of 
p, p, q and 5. In Table I we present some of our estimates 
for the exponent v at p — 0.75 and using in the simula- 
tions the conjectured value b c (S, x) giving in @- Our 
results are clearly consistent with the exponent v = ^ 
for all values of 5 and x- 

Up to now our results indicate that for arbitrary range 
of exclusion (size 5 = 0, 1, 2, . . . ) we can define appropri- 
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ate sliding-tag correlation functions exhibiting anoma- 
lous behavior. The existence of such anomalous behavior 
is connected to the fact that the dynamical behavior of 
these exclusion models are in the KPZ universality class. 
In fact the Bethe-ansatz solution |l6| of the quantum 
Hamiltonian associated to the asymmetric diffusion of 
particles with arbitrary size S (0,1,...) gives a dynami- 
cal critical exponent z = | in the KPZ universality class. 

A key ingredient in all the previous diffusion processes 
is that independently of the particle's sizes always their 
order are fixed (even for size particles). Now we are 
going to consider more general models where this order- 
ing requirement is partially brocken. These models will 
correspond to a generalization of the diffusion dynamics 
of size S particles, for the situation where the particles 
have a negative "size" S (S = —1,-2,...). In those 
models the particles are ordered from left to right and 
occupy as before the positions {y(n, t), n = 1, . . . , N p }), 
with the constraint y(n + l,i) > y(n,t) + S. Then for 
negative values of S the nth particle not necessarily is 
on the right of the (n — l)th particle, but can be at a 
distance —S on the left of this particle. This negative 
size S gives a measure of the ordering of the particles. 
In the limit S —>■ ~oo we have no ordering at all and the 
particles should behave as noninteracting random walks, 
with the absense of any anomalous behavior. It is inter- 
esting to note that the related interface model, given by 
the surface with cartesian coordinates {n,y(n,t)}, can 
exhibt hills of arbitrary height but valleys of maximun 
depth — S, in units of lattice spacing. 

In Fig. 3 we show some examples of our simulations for 
S = —1, p = 0.75,q = 0.25 and density p = 1 in a lattice 
of N s = 150 000 sites. As we can see in this figure an 
anomalous behavior is also detected at the critical value 
b c = 0.1667. This value indicate that the conjecture (||), 
that gives b c = | for the present case, is also valid for 
negative values of S. Similar results are also obtained for 
other negative values of S (see Table I) . 

Inspired by those last results we were able to extend 
the Bethe-ansatz calculations, already known |ll| for par- 
ticles with size S > 0, for the equivalent Hamiltonian de- 
scribing the particles with negative values of S jl7|. In 
the extreme anisotropic limit p = 1, q = the eigenener- 
gies of the associated Hamiltonian is given by 

E = J^(l- Zj )/2, (6) 

3=1 

where the roots {zj} are given by the Bethe-ansatz equa- 
tions 

(1 + ^"^1 ~ = ~2 N 'e-^ J] 

(7) 



where j = 1, . . . , N p , m = 0, 1, ... , N p — 1. The above 
equations are more difficult to handle analitically than 
the corresponding ones in the cases where S > 0, even 

TABLE I: Some values of the estimates of the exponent v 
of the correlation function (^) evaluated at p — 0.75, and for 
molecules of size S. These results were obtained at the tag 
velocity vt ag = b c pv p , with b c giving by (|) 







8 = -2 


S = -1 


S = 


S = 2 


X = 
X = 


1 

f 

5 


0.32 
0.31 


0.31 
0.32 


0.33 
0.31 


0.32 
0.33 



TABLE II: Finite-size estimates ^jv^.iYi of the critical dy- 
namical exponent evaluated by solving (ffl) for N s = N p . The 
extrapolated results are also shown. 



N S ,N' S 


50,52 


100,102 


150,152 


200,202 


CO 


S = -l 
S = -2 


1.61344 
1.72222 


1.56046 
1.62689 


1.54122 
1.58882 


1.53097 
1.56769 


1.500±0.001 
1.504±0.005 



at the " half-filling" x = \i where a simplification occurs 
@, § §. The real part of the energy gap G Ns = 
E\ — Eq behaves for large N s as Re(GN s ) ~ N~ z , pro- 
ducing for a pair of lattice sizes Ns, N' s the finite-size es- 
timate z Ns , N > = -lnRe(G Ns )/Re(G N >)/ln(N s /N' s ) for 
the dynamical critical exponent z. We solved numeri- 
cally (0) for S = -1 and S = -2 at the filling N p = N s . 
The estimates zn s ,n s +2 for some values of N s are pre- 
sented in Table II. We see clearly that the exponent tend 
toward the KPZ value z = |. 

In conclusion, our results indicate that as long some 
type of relative ordering exists (|5| < oo), always we 
are going to find anomalous behavior for the sliding-tag 
correlation (||) at the critical tag velocity vtag = pv p b c 
obtained from ([!]) and (|5|). This fact can also be under- 
tood in terms of the related interface growth model, and 
it implies that as long \S\ is finite, we should expect the 
long time fluctuations being governed by a KPZ model. 
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Figure Captions 



Figure 1 - The time dependence of the sliding-tag cor- 
relations (||) with size S — 2 particles for some values 
of b. The parameters of the Monte Carlo simulations 
are N s = 150 000, N p = 30 000, p = 0.75. The av- 
eraged estimated slope for the curves are 1.00 (6 = 1), 
0.90 (6 = 0.6), 0.84 (b = 0.3) and at b = b c we have the 
estimate 2v = 0.65 ± 0.03. 

Figure 2 - Same as figure 1 for size S = particles. 
The parameters of the Monte Carlo simulations are N s — 
120 000, N p = 30 000, p = 0.75. The averaged estimated 
slope for the curves are 0.97 {b = 1), 0.95 {b = 0.5), 
0.93 (b — 0.4) and at b = b c we have the estimate 2v — 
0.63 ±0.03. 

Figure 3 - Same as figure 1 for size S = — 1 parti- 
cles. The parameters of the Monte Carlo simulations 
are N s = 150 000, N p = 150 000, p = 0.75. The aver- 
aged estimated slope for the curves are 0.86 (b = 0.4), 
0.81 (b = 0.3), 0.76 (b = 0.1) and at b = b c we have the 
estimate 2v = 0.63 ± 0.03. 
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